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INTRODUCTION 


We  have  employed  BENCHMARK-76,  a version  of  the  AIRCHEM  code  with 
a 64  species  set  and  a nominal  496  reaction  set  to  obtain  species  den- 
sities as  a function  of  time  and  altitud|^for  |he  following  disturbed  con- 
ditions. The  prompt  ionization,  N = 10  cm  and  the  delayed  ioni- 
zation is  given  by 


Q(t)  = Qq  (1  + (1) 

8 3 

where  t is  the  time  in  seconds  and  = 10  ion-pairs/cm  -s.  The  rela- 
tive distribution  of  the  species;  N_^,  0-^,  N^,  0^,  N('^S),  N(^D),  0(^P) , 
11  1 z z 

0(  D) , 02(a  A)  and  02(b  Z)  produced  by  ionization  under  disturbed  con- 
ditions is  taken  to  be:  0.64,  0.16,  0.14,  0.06,  0.45,  0.61,  1.28,  0.10, 
0.24  and  0.05  particles  per  ion  ga^r,  respectively.  This  distribution 
closely  follows  that  of  Gilmore,  ’ except  that  [02(^A)]  is  taken 
as  0.24  and  [0]  is  1.28.  Cases  for  daytime  and  nighttime  conditions  are 
reported;  the  calculations  correspond  to  about  a three  hour  interval 
centered  around  noon  and^midnight.  Appropriate  neutral  densities  have 
been  outlined  elsewhere.  The  program  coding  and  fixed-input  parameters 
have  been  previously  discussed.  Results  were  obtained  from  the  BENCH- 
MARK-76 program  for  the  altitude  region  10-80  km  inclusive  for  daytime 
conditions  and  30-80  km  for  nighttime. 

It  is  the  purpose  of  this  report  to;  (1)  display  selected  computed 
results  as  a function  of  time  and/or  altitude,  (2)  show  the  results  of 
comparison  and  sensitivity  studies  of  a very  limited  nature,  (3)  report 
the  temporal  and  altitude  variations  of  the  computed  equivalent  rate 
coefficients  (ERC's),  and  (4)  indicate  the  type  of  detailed  analysis 
available  from  the  BENCHMARK-76  output.  ( A copy  of  these  BENCHMARK-76 
results  may  be  obtained  by  writing  the  authors.) 


^E.  L.  Lortie,  M.  D.  Kregel  and  F.  E.  Niles,  "AIRCHEM:  A Computational 
Technique  for  Modeling  the  Chemistry  of  the  Atmosphere,"  BRL  Report  No. 
1913,  August,  1976.  (AD#  A030157) 

2 

F.  Gilmore  as  quoted  by  B.  F.  Myers  and  M.  R.  Schoonover,  "Electron 
Energy  Degradation  in  the  Atmosphere:  Consequent  Species  and  Energy 
Densities,  Electron  Flux,  and  Radiation  Spectra,"  DNA  3513T,  3 Jan 
75,  Table  6. 

^J.  M.  Heimerl  and  F.  E.  Niles,  "BENCHMARK- 76 : Model  Computations  for 
Disturbed  Atmospheric  Conditions  I.  Input  Parameters,"  BRL  Report 
No.  2022,  October  1977.  (AD  #A050355) 

4 

F.  E.  Niles  and  J.  M.  Heimerl,  "Selected  Neutral  Species  Profiles 
0-100  km,"  BRL  Memo  Report  No.  2767,  July  1977.  (AD#  A042620) 
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RESULTS 


A.  Selected  Densities 


The  electron  density  computed  as  a function  of  time  is  shown  for 
each  decade  of  altitude  for  day  and  night  conditions  in  Figures  1 and  2, 
respectively.  Q,  the  electron  production  as  a function  of  time  is  shown 
at  the  uppermost  section  of  each  figure.  The  smooth  curves  were  drawn 
between  computed  values  at  1,  2,  4 and  7 units  of  each  decade  in  time. 

Except  at  80  |m  during  the  night  and  above  60  km  during  the  day  it 
appears  that  by  10  seconds  the  electron  decay  follows  the  production 
function  decay.  Three-body  attachment  is  the  dominant  electron  depopu- 
lating process  at  late  times  (except  at  70  km  and  80  km)  and  Q is  the 
dominant  production  mechanism  (except  for  late  times  at  80  km  during  the 
day) . The  reason  for  this  peculiarity  is  discussed  below  (Section  E)  as 
an  example  of  the  detailed  analysis  available  from  the  BENCHMARK-76  output. 

Some  of  the  curves  cease  abruptly  in  these  (and  other)  figures 
because  computations  are  ai^tomatically  halted  whenever  the  electron 
density  falls  below  10  cm  , the  assumed  quiescent  background  electron 
density. 

Figures  3 and  4 show  the  sum  of  the  negative  ions  (exclusive  of 
electrons)  computed  as  a function  of  time.  Conditions  correspond  to  those 
in  Figures  1 and  2,  respectively.  Figures  5 and  6 show  the  corresponding 
plots  for  the  sum  of  the  positive  ions.  This  sum  is  insensitive  to 
day/night  conditions.  To  avoid  crowding,  odd  decades  of  altitude  have 
been  plotted  in  Figure  5 and  even  decades  in  Figure  6. 

B.  Comparison 

Figures  7 and  8 show  the  electron  density  computed  as  a function  of 
time  at  60  km  altitude  fgr  both  the  BENCHMARK-76  results  (BRL-76)  and  the 
DCIIEM  results  of  Scheibe  for  daytime  and  nighttime,  respectively.  Figures 
9 and  10  show  similar  comparisons  for  the  total  positive  and  negative 
ions  under  daytime  and  nighttime  conditions,  respectively.  The  agree- 
ment is  very  good  to  excellent.  These  two  co^es  employ  different  integra- 
tion schemes  and  different  overall  chemistry.  However,  a detailed  exami- 
nation of  the  controlling  chemistry  at  60  km  might  show  strong  similarities. 
Such  an  examination  is  beyond  the  ^cope  of  the  present  report.  Comparisons 
made  for  other  values  of  Q and  N show  basic  agreement  of  results  of 
these  two  codes.  We  take  ?hese  comparisons  to  indicate  that  the  agree- 
ment is  real  rather  than  fortuitous. 

^W.  S.  Knapp,  "A  Simplified  D-Region  Chemistry  Model  for  Nuclear  Environ- 
ments," DNA-2850T,  April  72;  also  M.  Scheibe,  Mission  Research  Corpora- 
tion, private  communication,  1976. 

0 

For  example.  The  BENCHMARK-76  reaction  set  includes  neutral  odd-hydrogen 
reactions,  while  Scheibe' s set  does  not. 

.1.  M.  Heimerl  and  F.  E.  Niles,  "Br.NCHMARK-76 : Model  Computations  for 
Disturbed  Atmospheric  Conditions  III.  Results  for  Selected  Excitation 
Parameters  at  60  km,"  BRL  Technical  Report  ARBRL-TR-02051 , March  1978. 
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Figure  1.  Logarithm  of  the  computed  daytime  electron  density  vs. 

logarithm  of  time.  Parameters  are  the  altitudes  in  km. 

Top  curve  is  the  ion-pair  production  as  a function  of 
time.  ,Q  is  the  delayed  ionization  parameter,  in  ion  pairs 
cm"^s‘'.°  N is  prompt  ionization  parameter  (electron 
density)  in°cm"^. 
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Figure  2.  Same  as  Figure  1,  except  nighttime  conditions. 
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Figure  3.  Logarithm  of  the  sum  of  the  computed  negative  ion  densities 
vs.  logarithm  of  time  for  daytime  conditions.  Parameters 
are  the  altitudes  in  km. 
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Figure  4,  Same  as  Figure  3,  except  nighttime  conditions. 
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Figure  5.  Logarithm  of  the  sum  of  the  computed  daytime  positive  ion 
densities  vs.  logarithm  of  time.  Parameters  are  the  odd 
altitudes  in  km.  (See  text.) 

15 


POS  lONSj 


Figure  6.  Same  as  Figure  5,  except  nighttime  conditions  for  the  even 
altitudes  (See  text). 
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Figure  7.  Jfe  daytime  electron  density,  at  60  km.  computed  by  BRL's 
BENCHMARK-76  code  (•)  and  Scheibe's  DAIRCHEM  code  (x). 
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Figure  8.  Same  as  Figure  7 except  nighttime  conditions. 
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Figure  9.  The  total  positive  ion  densities  (— ) and  total  negative 
ion  densities  { — ) at  60  km  for  daytime  conditions  as 
computed  by  BRL's  BENCHMARK-76  code  (•)  and  Scheibe's 
DAIRCHEM  code  (x).  ,q 


cm 


Figure  10.  Same  as  Figure  9,  except  nighttime  conditions. 
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C.  Sensitivity 


Of  the  210  ion-ion  recombination  reactions  used  in  the  BRL  code 
only  8 have  been  measured  in  ^he  labora^ojy.  jThi  remaining  202  were 
assigned  a value  of  2.0  x 10  (T/300)  cm  s . Since  there  is  some 

discussiongOven  as  to  the  order  of  magnitude  of  some  of  the  measured 
reactions,  computations  were  made  for  the  two  cases  wherein  all^  of^the 
ion-ion  r^ipombination  ratp  coefficients,  n.'s,  are  multiplied  by  10“ 
and  by  10  . The  results  of  this  exercise,  for  nighttime  conditions  at 

an  altitude  of  60  km  are  shown  in  Figures  11,  12  and  13  for  electrons 
negative  ions,  and  positive  ions,  respectively.  Though  there  are  large 
changes  in  both  the  total  positive  ion  and  total  negative  ion  densities, 
the  change  in  the  electron  density  barely  exceeds  a factor  of  two  for  a 
factor  of  100  change  in  a^.  This  insensitivity  to  the  precise  value  of 
rate  coefficients  can  be  expected  whenever  the  production^gf  electrons 
is  governed  by  the  source  term  rather  than  the  chemistry.  For  electrons 
it  appears  as  though  it  is  of  the  utmost  importance  to  include  the 
correct  processes,  rather  than  have  very  accurate  rate  coefficients.  On 
the  other  hand,  the  hundredfold  range  of  values  assumed  for  the  unmeasured 
ion-ion  recombination  coefficients  can  make  substantial  difference  in  the 
range  of  values  computed  for  both  the  total  negative  ion  and  the  total 
positive  ion  densities,  as  can  be  seen  in  Figures  12  and  13  respectively. 

We  shall  now  show  that  a tenfold  spread  in  each  of  these  values  is  expected 

Whenever  [N  ] >>  [e] 

q - «gff[N'][N^]  , (2) 

where  q is  the  production  of  negative  ions,  is  the  effective  ion-ion 

3 — 

D.  Smith  and  M.  J.  Church,  "Binary  Ion-Ion  Recombination  Coefficients 
Determined  in  a Flowing  Afterglow  Plasma,"  Int.  J.  Mass  Spec,  and 
Ion.  Phys.  19,  185,  1976. 

9 

For  the  case  where  selected  a.'s  are  varied  by  the  factors  see  "The 
Effects  of  Uncertainties  in  tfte  Two-Body  Ion-Ion  Recombination  Coeffi- 
cient upon  Computed  Ion  Distributions  in  the  Stratosphere  and  Mesosphere,' 
by  J.  M.  Heimerl,  BRL  Report  No.  1972,  March  1977  (AD#  A038263) . Table 
2 of  this  reference  contains  three  typographical  errors  and  one  com- 
putational error.  At  60  km  the  fourth  reaction  should  read  H (H-O).  + 

NOj  ; at  40  km  the  third  reaction  H (H20)^  + £0j”(H20)  and  at  30  km"^ 

the  first  reaction  H^(H20)^  + CO^  (H2O) . The  values  in  Table  3 of  this 
reference  then  correspond  to  variations  in  the  corrected  ion-ion 
recombination  reactions  listed  in  Table  2^  The  computations  at 
40  km  should  have  been  carried  out  with  H (H20)^  + N02"(H20)  as  the 
second  reaction.  Recomputation  shows  some  changes  in  the  percentages 
in  Table  3 but  the  conclusions  therein  remain  unchanged. 

^^For  a counter  example  wherein  chemistry  dominates  the  production  terra, 
see  F.  E.  Niles,  "Electron  Density  in  the  Beta  Patch  Following  a Spartan 
Burst,"  BRL  Report  No.  1518,  December  1970,  Figure  9.  (AD#  879464) 
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Figure  11.  Sensitivity  of  the  computed  nighttime  electron  density  at 
60  km  to  simultaneous  decade  variations  in  each  value  of 
the  210  ion-ion  recombination  coefficients  used  in  the 
BENCHMARK-76  code. 
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Figure  12. 


Sensitivity  of  the  computed  nighttime  total  negative  ion 
density  for  the  conditions  stated  in  Figure  11. 
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Figure  13.  Sensitivity  of  the  computed  nighttime  total  positive  ion 
density  for  the  conditions  stated  in  Figure  11. 
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recombination  coefficient  and  [N  ] and  [N^]  are  the  total  negative  and 
total  positive  ion  densities,  respectively,  Undej;  steady-state  condi- 
tions, i.e,,  d[N"]/dt  = 0,  and  whenever  [N"]  =“  [N  ] , equation  (2) 
reduces  to 

q or  (3a) 

q^ctefft^’l^  • 

For  the  results  shown  in  Figures  12  and  13  the  production  of  negative 
ions  remains  unaltered  and  so  q = constant.  Since  a ® Ea.  a factor  of 

100  change  in  every  is  reflected  as  a factor  of  100  change  in  ag££- 

For  such  a change  equations  (3a)  and  (3b)  show  that  both  [N  ] and  [N~] 
will  change  by  a factor  of  ten. 

This  factor  of  ten  is  observed  at  the  latest  times  in  Figures  12 
and  13  because  the  conditions  under  which  this  factor  applies,  i.e,, 

[N  ] >>  [e] , [N  ] =»  [N^]  and  steady  state,  are  all  found  when  t > lO^s. 
(Compare  Figures  11  and  12  and  Figures  12  and  13.  In  Section  D we  shall 
find  that  the  steady  state  condition  is  reached  by  about  one  second.) 

D.  Equivalent  Rate  Coefficients 

Once  the  individual  species  concentrations  have  been  computed  as  a 
function  of  time  and  altitude  a set  of  equivalent  rate  coefficients, 

ERC ' s , can  be  generated . 

The  effective  electron-ion  recombination  coefficient  is  defined  by 

«d  = ^Rei/[e]J:[N^^]  (4) 

where  SR  ^ is  the  sum  of  the^rates  for  electron-ion  recombination,  [e] 
the  electron  density  and  [Nj^  ] is  the  density  of  the  k positive  ion 
species.  This  sum  is  over  all  positive  ion  species. 

There  are  24  reactions^ ^ that  make  up  the  ER  . of  Equation  (4).  A 
suffijienj  number ^of  reactions  is  carried  in  the  model  so  that,  except 
for  N , 0 and  O2  (N2) , each  positive  ion  species  jeacts  wit^  an  electron. 
(The  two  products  of  the  NO  + e reaction,  viz,  N(  S)  and  N(  D) , are  carried 
explicitly  as  two  separate  reactions.) 

Figures  14  and  15  show  the  temporal  variation  of  for  daytime  and 
nighttime  conditions,  respectively.  The  numbers  on  the  curves  refer  to 
altitude,  z,  in  kilometers.  For  the  lower  altitudes  approaches  limiting 
values  determined  by  the  relative  composition  of  cluster  ions  and  their 
respective  electron-ion  recombination  coefficients.  At  the  higher 
altitudes  the  limiting  values  of  are  not  yet  attained,  even  after  10^ 
model  seconds. 

^^Reactions  138-162  inclusive;  see  Appendix  for  reaction  list. 
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log  Ms) 

Figure  14.  Effective  electron-ion  recombination  rate  coefficient 
as  a function  of  time  for  daytime  conditions.  Parame 
are  the  altitudes  in  km. 


In  figure  14  the  crossover  of  a.  for  50  and  60  km  is  explained  by 
differing  cluster  composition  with  altitude^  To  be  specific,  consider 
tlje  value  of  a,  at  50  and  60  km,  both  at  10  model  seconds.  At  50  km 
accounts  for  64%  of  the  total  positive  ions,  NO  (H2O)  - 20%, 

NO^CH^O)^  - 9%  and  - 8%.  The  corresponding  electron  ion-recombi- 


nation coefficients,  . employed  in  the  code  are:  5.20,  1.02,  2.04 
and  6.22  x 10  ^ cm^/s,  respectively.  Noting  that  R . = a • [N, [e] , we 
find  from  Equation  (4)  that  a^(z=50  km,  t=10  s)  = 4.2  cm'^/s.  At  60  km, 
accounts  for  40%  of  the  total  positive  ions,  - 26%, 

N0^(H,0)^  - 24%  and  NO^(H-O)  - 9%.  The  corresponding  a . are:  5.29, 

— *”  - o 3 ^ ^ 

6.33,  2.08  and  1.04  x 10  cm  /s.  Again  from  Equation  (4)  we  find  that 

4 3 -4 

a^tz=60  km,  t=10  s)  = 4.4  cm  /s.  In  a similar  fashion  the  "dip"  at  10 


seconds  for  30  km  can  be  explained  by  differing  positive-ion  composition 
with  time. 


The  effective  electron  attachment  rate.  A,  is  given  by 

A = ER^/[e]  (5) 

12 

where  ZR^  is  the  sum  of  the  rates  for  attachment  processes  and  [e]  is 
the  electron  density.  Figures  16  and  17  show  values  of  A plotted  against 
time  for  daytime  and  nighttime  conditions,  respectively.  The  numbers  to 
the  right  of  the  curves  denote  the  altitude,  z,  iy  kilometers.  The 
variations  observed  in  the  values  of  A for  t > 10  s and  for  z > 60  km 
will  now  be  shown  to  reflect  variations  in  the  ozone  density. 


ERa 


From  an  examination 
*^124  '^125  *^126 


O2  ^ 


and 


°2  ^ 


of  the 
where 

output 

it  is 

found  that 

O3  ^ e 

->• 

o" 

^ °2 

R124* 

e ^ O2 

->• 

0" 

^ °2 

R125 

e + N^ 

O2 

- ^ N,. 

R126 

Substituting  this  approximation  into  equation  (5)  we  find  that 

A = C(1  + B)  (6) 


where 


^ " ’^125^°2^  ^126^^2^  f^2^ 


(7) 


B = kj24[0^]/c  . (8) 

_ 

Keactions  122-126  inclusive;  sec  Appendix  for  reaction  list, 

* 

Read  as  Reaction  124.  A complete  listing  of  reactions  and  references  may 
be  found  in  the  Appendix. 
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For  a given  altitude,  C defined  by  equation  (7)  is  approximately  constant. 
Since,  again  at  a given  altitude,  the  value  of  ^,2^  is  fixed,  B in  equa- 
tion (8)  reflects  the  changes  in  [0^] . This  reflection  shows  up  in  the 
expression  for  A,  equation  (6) . 

The  value  of  B in  equation  (8)  defines  the  importance  of  reaction 
124  relative  to  the  three-body  attachment  reactions,  i.e. 

^ ~ ^124^~^125  * *^126^  ■ i ^ show  the  values  of  B for  daytime 

and  nighttime  conditions,  respectively.  From  Table  2 we  find  that  B is 
small  or  constant  in  time,  or  both  (1)  for  altitudes  of  50  km  or  less 
over  the  entire  time  interval  studied  or  (2)  for  the  altitude  region  60- 
80  km  and  for  t < lO^s.  From  equation  (6)  we  expect  the  value  of  A to  be 
constant  in  time  for  these  conditions  and  this  js  observed  in  Figure  17, 

For  the  altitude  region  60-80  km  and  for  t > 10  s the  values  of  B in 
Table  2 are  seen  to  vary  with  time.  The  corresponding  variations  in  A 
are  again  found  in  Figure  17. 

A similar  analysis  follows  for  the  values  of  B for  daytime  condi- 
tions given  in  Table  1. 

We  will  now  briefly  and  qualitatively  examine  the  proximate  cause 
of  the  [Oj]  variations.  An  examination  of  the  output  reveals  that  for 
the  higher  altitudes  by  t = 10~^s  the  production  of  [0^]  is  almost  entirely 
due  to  the  reaction 


0 + 02+M->-02+M.  R117 

This  reaction  shows  that  at  a given  altitude,  i.e.  both  [O2]  and  [M]  are 
constants,  the  ozone  production  will  track  the^atomic  oxygen  density.  The 
[0]  is  found  to  be  constant  for  times  up  to  lO'^s  at  the  higher  altitudes; 
thus  [Oj]  production  is  approximately  constant.  The  ozone  density  can 
rise,  however,  ^ecause  ^he  net  loss  rate  of  ozone  due  to  reactions  with: 

H,  N,  NO,  0,  0(  D) , and  electrons  (corresponding  to  reactions  70, 

81,  88,  93,  94,  100,  104  and  124,  respectively)  falls.  These  higher 
[Oj]  values  are  not  maintained,  because  by  lO^s  the  atomic  oxygen  density 
declines  leading  to  the  eventual  decline  in  [Oj] . 

The  effective  electron  detachment  rate  is  defined  by 

D = 2:Rj^/Z[N^‘].  (9) 

where  th|j^Rp  are  the  rates  for  individual  detachment  processes  and  [N. "] 
is  the  k^  negative  ion  density.  The  sums  are  over  all  detachment  rates 
and  all  negative  ion  densities,  respectively.  The  two  detachment  processes 
considered  in  the  code  are  photodetachment  and  two-body  collisiom^l  detach- 
ment; the  former  is  a set  of  six  reactions  and  the  later  eleven. 


13 

Reactions  1-5  and  11  for  photodetachment;  127-137  inclusive  for  two 
body  collisional  detachment;  see  Appendix  for  reaction  list. 
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TABLE  1.  VALUES  OF  B FOR  DAYTIME  CONDITIONS 


Figures  18  and  19  show  D vs  time  for  daytime  and  nighttime  condi- 
tions, respectively.  Altitudes,  in  kilometers,  are  listed  nearby  their 
corresponding  curves.  In  both  figures  the  most  notable  feature  is  that 
no  curve  attains  ^ constant  value  at  late  times.  ("Late  time"  for  the 
10  km  case  is  10  s.)  A detachment  process  that  is  important  at  all 

the  altitudes  studied  is 


0"  + N2  N2O  + e.  R129 

This  reaction  has  been  assigned  the  upper  limit  of  1 x 10  cm^/s.^^ 

To  examine  the  effect  this  value  might  have  on^^Jhe  |ffe^|ive  electron 
detachment  rate  we  assigned  a value  of  1 x 10  cm  /s  and  exercised 
the  model  over  the  altitude  range  50-80  km  for  daytime  conditions.  Table 
3 shows  the  comparison  of  the  values  of  D as  a function  of  time.  Since 
R,29  “ ^129’  equation  (9)  predicts  that  values  of  D should  be  greater 
the  larger  the  value  of  k.-^.  This  behavior  is  observed  in  Table  3. 

As  may  be  seen  in  this  table  the  percentage  differences  in  the  value  of 
D in  some  cases  amounts  to  about  50%  for  a hundredfold  change  in  the  value 
of  kj29’  Another  result  of  this  exercise  was  that  while  reaction  (129) 
was  no^longer  dominant,  it  was  still  significant  and  could  not  in  general 
be  ignored. 

The  effective  recombination  coefficient,  i|;,  is  defined  by 

= (1  X)  (Uj  + Xa^) , (10) 

where  X is  the  ratio  of  the  total  negative  ion  density  to  the  electron 
density  (other  terms  have  been  define^^.  In  the  steady  state,  a criterion 
many  times  employed  in  hybrid  models,  should  equ^l  the  source  term 
divided  by  the  square  of  the  electron  density,  Q/[e]  . At  a gi^en  altitude, 
the  value  of  ip  (solid  line  in  Figure  20), and  the  value  of  Q/[e]  (dashed 
line  in  Figure  20)  are  about  equal  (i.e.,  within  about  20%)  for  times 
greater  than  about  one  second  (see  arrow) . This  result  also  holds  true 
for  daytime  conditions  (not  shown). 

E.  Example  of  Detailed  Analysis 

In  this  section  we  shall  illustrate  the  type  of  detailed  analysis 
available  from  the  BENCHMARK-76  program  output.  In  Figure  1 the  electron 
density  r*  80  km  does  not  follow  the  source  term,  Q,  at  late  times,  but 
actually  uegins  to  increase.  This  is  explained,  as  we  shall  see,  by  an 
additional  source  for  electrons. 

T4 

W.  Lindinger,  D.  L.  Albritton,  F.  C.  Fehsenfeld  and  E.  E.  Ferguson, 
"Reactions  of  0 with  N2,  N2O,  SO2,  NH^,  CH.  and  and  C2H2  with 

0 from  300°K  to  relative  kinetic  energies  or  ~ 2 eV,"  J.  Chem.  Phys. 

W,  3238-3242,  1975. 

^^F.  C.  Fehsenfeld,  private  communication,  1977. 

^^See  for  example,  A.  P.  Mitra,  "D-Region  in  Disturbed  Condition, 

Including  Flares  and  Energetic  Particles,"  J.  Atmos,  and  Terr.  Phys. 

37,  895-913,  1975. 
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TABLE  3.  COMPUTED  VALUES  OF  D (s‘^)  AS  A FUNCTION  OF  TIME  AND  VALUE  OF 
FOR  FOUR  ALTITUDES. 


ALT 

*^12^ 

(cc/s) 

time 

\^(s) 

10^ 

10^ 

80 

10-'^ 

I.32(+2)* 

I.39(+2) 

1.39(+2) 

9.16(+I) 

5.0I(+1) 

80 

10-^' 

9.57(+l) 

9.53(+l) 

8.89(+l) 

6.56(+l) 

3.85(+I) 

70 

10-^' 

1.4I(+2) 

1.44(+2) 

I.48(+2) 

9.67(+l) 

3.82(0) 

70 

10-'^ 

9.46(+l) 

8.99(+I) 

7.38(+I) 

5.57(+I) 

3.12(0) 

60 

10-^' 

I.Il(+2) 

8.43(+l) 

I.94(+I) 

1.92(0) 

I.02(-l) 

60 

10-^4 

7.34(+I) 

5.25(+I) 

I.I7(+I) 

1.69(0) 

9.23(-2) 

50 

10-'2 

2.66(+I) 

5.90(0) 

6.06(-I) 

I.43(-l) 

I.34(-2) 

50 

10-^4 

1.44(+1) 

3.29(0) 

4.24(-l) 

9.61(-2) 

5.56(-3) 

+2 


Read  as  1 . 32  x 10 
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Figure  20.  Composite  plot  of  the  logarithm  of  the  total  recombination 
coefficient,  ijj  (solid  line),  and  of  Q/[e]^  (dashed  line) 
as  a function  of  the  logarithm  of  time  for  nighttime 
conditions.  Parameters  are  the  altitudes  in  km.  The 
arrow  points  to  decade  after  which  =»  Q/[e‘^]. 


Figure  21  shows  a log-log  plot  of  the  electron  populating  rates  for 
specific  reactions  as  a function  of  time.  These  reactions  are  selected 
from  a host  of  electron  populating  reactions  listed  in  the  BENCHMARK-76 
output,  on  the  basis  that  at  least  at  some  decade  in  time,  the  specific 
rate  is  > 10%  of  the  maximum  rate  at  that  decade. 

Over  most  of  the  time  coordinate^electrons  are  populated  by  the 
source  term,  Q.  But  at  times  near  10  the  reaction 

NO  + hv  -»■  NO^  + e R51 

becomes  the  dominant  source  of  electrons.  Figure  22  shows  the  corres- 
ponding plot  for  electron  depopulation.  This  plot  shows  that  electron- 
ion  recombination  is  the  dominant  electron  loss  process  during  all  model 
times  for  this  altitude. 


SUMMARY 

BENCHMARK-76  has  been  employed  to  obtain  species  densities  as  a 
function  of  time  and  altitude  for  the  followj^g  excitation  conditions. 

The  prompt  ionization,  N^,  is  taken  t^  ^e  10  cm  and  the  delayed 
iogization  is  giv^n  b^  Q = Qj(l  t)”  ’ . where  t is  in  seconds  and  = 

10  ion  pairs  cm  -s  . Calculations  were  made  for  three  hour  intervals 
around  noon  (DAY)  over  the  altitude  range  10-80  km  and  midnight  (NIGHT) 
over  the  altitude  range  30-80  km.  Analysis  of  these  results  is  summarized 
as  follows. 

First,  at  sufficiently  late  times  and  at  all  but  the  highest  alti- 
tudes, the  computed  electron  density  tracks  the  delayed  ionization.  The 
computed  total  negative  ion  density  is  a strong  function  of  altitude,  as 
expected.  The  computations  show  that  both  the  total  positive  ion  and 
the  total  negative  ion  densities  are  insensitive  to  DAY/NIGHT  conditions. 

Second,  the  comparison  between  the  BENCHMARK-76  (BRL)  and  UCHEM  (MRC) 
computations  of  the  electron  density,  the  total  negative  ion  and  the 
total  positive  ion  densities  at  60  km  is  good  to  excellent. 

Third,  the  overall  effect  of  simultaneously  varying  every  ion-ion 
recombination  coefficient  a hundredfold  produced  a maximum  twofold 
change  in  the  electron  density  and  the  expected  tenfold  change  in  both 
total  positive  and  total  negative  ion  densities. 

Fourth,  the  effective  electron-ion  recombination  coefficient,  a^, 
is  found  to  be  a strong  function  of  altitude  and  time.  For  altitudes 
greater  than  40  km  i|  is  found  to  be  still  approaching  its  quiescent 
value,  even  after  10^  model  seconds. 
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DAY 


log  f (s) 


Figure  21.  Logarithm  of  electron  populating  rates  at  80  km  vs.  the 
logarithm  of  time  for  daytime  conditions.  Specific 
reactions  are  identified,  Q is  the  delayed  production 
of  electrons. 
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Figure  22.  Logarithm  of  electron  depopulating  rates  at  80  km  vs. 

logarithm  of  time  for  daytime  conditions.  Specific 
reactions  are  identified. 
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Fifth,  the  effective  electron  attachment.  A,  is  dominated  by  three 
body  attachment  to  0^  and  is  not  a function  of  time  for  altitudes  below 
about  60  km.  Increases  of  up  to  a factor  of  three  in  A at  late  times 
and  at  altitudes  higher  than  60  km  are  found  to  be  caused  by  an  increase 
in  ozone  which  permits  the  reaction  0^  + e ->■  0 + O2  to  compete  with  O2 

attachment. 

Sixth,  the  effective  electron  detachment  rate  is  found  to  be  a 
strong  function  of  time  and  altitude  and  under  no  conditions  employed 
did  it  reach  a constant  value.  The  effective  electron  detachment  rate 
is  dominat^^  by^the  reaction  0 + N2  ^ N^O  + e when  a rate  coefficient 
of  1 X 10  cm  /s  is  used.  Even  when  this  rate  coefficient  is  reduced 

a hundredfold,  the  reaction  in  general  can  not  be  neglected. 

Finally^  the  comparison  of  the  effective  recombination  coefficient, 
4),  and  Q/[e]  showed  that  a steady  state  is  reached  at  all  altitudes 
studied  at  about  one  second. 
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appendix 


The  table  in  this  Appendix  displays  the  nominal  496  reactions  that 
link  the  64  species  within  the  BENCHMARK-76  code.  Reactions  proceed  from 
left  to  right  only.  Each  rate  coefficient,  k,  is  constructed  from  the 
reaction  parameters  A,  B and  C,  i.e. 

k = A(T/300)®  exp  (-  C/T)  , 

where  T is  the  temperature  in  degrees  Kelvin.  Except  for  the  first  60 
reactions,  involving  a photon  as  a reactant,  the  reference  for  the  value 
of  the  rate  coefficient  is  given  to  the  right  of  the  parameter  C column. 
Further  details  can  be  obtained  from  Reference  3. 
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TABLE  A.  BENCHMARK-76  REACTION  SET 
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